Understanding Kondo Peak Splitting and the Mechanism of Coherent Transport in a 
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The peculiar behavior of Kondo peak splitting under a magnetic field and bias can be explained 
by calculating the nonequilibrium retarded Green's function via the nonperturbative dynamical 
theory (NDT). In the NDT, the application of a lead-dot-lead system reveals that new resonant 
tunneling levels are activated near the Fermi level and the conventional Kondo peak at the Fermi 
level diminishes when a bias is applied. Magnetic field causes asymmetry in the spectral density and 
transforms the new resonant peak into a major peak whose behavior explains all the features of the 
nonequilibrium Kondo phenomenon. We also show the mechanism of coherent transport through 
the new resonant tunneling level. 
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After the observation of the equihbrium Kondo 
phenomenon in a single-electron transistor (SET)!!, 
nonequihbriuni Kondo phenomenon has rapidly evolved 
into one of the highly debated subjects in condensed 
matter physics. As the phenomenon entails two the- 
oretically challenging field of study, namely nonequi- 
librium and strong correlation, thus far no theoretical 
study has been successfully able to explain the experi- 
ments exhibiting the nonequilibrium Kondo phenomenon 
fully[2, 3]. Therefore, a theoretical understanding of 
this phenomenon would be an essential development in 
the advancement of condensed matter physics. The 
most attractive aspect of the nonequilibrium Kondo phe- 
nomenon is the splitting of the Kondo peak under a 
magnetic field 2, 3]. According to the experiment per- 
formed by Amasha et a/. [3], nonequilibrium Kondo phe- 
nomena can be summarized as follows: (i) splitting vs. 
magnetic field, which is expressed by a simple relation 
Ak,s = g + IqI^bB, where Ak,s denotes half of the 
gap between the split Kondo peaks, the superscript 
denotes field-independence, and the subscript S denotes 
particle-hole symmetric case; (ii) splitting vs. gate volt- 
age, which is given by = A^ g - Co(Vg - V"g,o)^, 
where V^^o denotes the gate voltage in the middle of 
the Coulomb valley and the curvature of the parabola 
Co is field- independent; (iii) splitting vs. Kondo tem- 
perature, which appears to be logarithmically decreas- 
ing, i.e., A^{Tk) oc — CiIuTk, where Ci is also field- 
independent; and (iv) the maximum Kondo temperature 
7km vs. the critical magnetic field Bq at the split- 
ting threshold, which is nonlinear. One notable issue 
common to all the abovementioned features is the field- 
independent behaviors of A^^, Co, and Ci. Existing 
theories [3| cannot explain the field-independent behav- 
iors in (i), (ii), and (iii) and the nonlinear behavior in 
(iv). We will show that these can be explained by the 
nonperturbative dynamical theory (NDT)0]. 

The result of NDT shows that the Kondo peak in equi- 
librium comprises two different types of coherence whose 
spectral weights will be denoted by Zg^ and Zg'^ below. 
In the first type of coherence, the spin-down electron. 
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FIG. 1: Motions of spin-down electron in the first (a) and 
second (b) type of coherence. In (b), left one corresponds to 
Z'^^ while right one to Z'^'^. 



for instance, travels back and forth through the dot as 
shown in Fig. 1 (a), while it moves as that shown in 
Fig. 1 (b) in the second type of coherence. The spin- 
up electron arrives toward the dot from the left or right 
lead in both the cases. The spectral density has a sin- 
gle resonant peak in equilibrium as shown in the inset of 
Fig. 2. Because of two separate metallic reservoirs, the 
maximum of TrApd] {ijj) is 0.5 instead of unity. 

The most interesting result of NDT appears when a 
bias is applied. A part of the spectral weight of Zg^ 
transfers from the Fermi level to the new resonant tunnel- 
ing levels near the Fermi level, while the spectral weight 




FIG. 2: The spectral density of the SET for U = lOA under 
a bias, where A = r/2. The Zeeman shift is 0.5A. The inset 
is for that under equilibrium. 
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Zg^, which corresponds to the second type of coherence, 
remains at the Fermi level. However, the weight Zg'" at 
the Fermi level is suppressed by decoherence due to bias. 
Interesting point is that the position of the new reso- 
nant level is independent of the applied field. This field- 
independence is responsible for the field independent fea- 
tures of the Kondo-peak splitting mentioned above. 

When magnetic field is applied, an asymmetry occurs 
in the spectral density in addition to the usual Zeeman 
shift IglnsB. As a result of asymmetry, one of the new 
resonant peaks becomes a major peak as shown in Fig. 
2. All the abovementioned features of the Kondo peak 
splitting phenomenon can be explained by the position 
of this major peak. In order to draw the spectral den- 
sity pdi (w) , we artificially choose the parameters of the 
spectral density to show the asymmetry and the suppres- 
sion of the central peak explicitly. The correct values of 
the parameters may be determined by the self-consistent 
calculation proposed in Ref. [5]. 

The result of NDT gives the positions of the new 



resonant peak as ±^Zg^U/2 in the Kondo regime 

with particle-hole symmetry, where U is the amount 
of Coulomb repulsion at the dot. Therefore, the field- 
independent splitting in feature (i) mentioned earlier is 



given by AjJ; — ^JZg^U/2. Since asymmetry is in- 
volved in the part that becomes U in the symmetric case, 
we propose the following expression of for the asym- 
metric case: 
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where Z = Z^^ + Z^^ and D = ZsF exp[7rC//8V2r]. 
Equation (1) recovers A^ ^ ii Z = Zs- If we use the 
notations given in Ref. [3], the wavefunction renormal- 
ization Z can be rewritten as Z = Z5exp[x(t4~^g,o)^] = 
ZsTk/Tk,Oj where Vg denotes the gate voltage and Tk.o 
is the Kondo temperature at gate voltage in the middle of 
the Coulomb valley Vg^. Then, the asymmetric behavior 
of A^ is given by the parabolic or logarithmic form 
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These expressions simultaneously explain all the above- 
mentioned features of the Kondo peak splitting in a qual- 
itative manner, except for the feature of the splitting 
threshold. The curvature of the parabola and the co- 
efficient in the — In TK-dependence are given by Cq = 
-Sv^FA^ gX/TrC/ and Ci = 8V2rA5^ ^/Trf/, respec- 
tively. Since the values of the constants Co and Ci are 
given by Co = -0.22/xeV/(mV)2 and Ci = ll^eV, re- 
spectively, using the experimental values F — 330/ieV, 
U = 1.2meV, x = 0.020(mV)-2, and AjJ, g = ll^eV, 
quantitative agreement with the experimental data is 
perfect. 



The threshold of the Kondo peak splitting depends on 
both the width of the peak and the positions of the major 
peaks produced by both the spin-up and spin-down elec- 
trons. Since the width of Kondo peak is linearly propor- 
tional to the Kondo temperature, the linear deviation of 
the Kondo temperature from its minimum, i.e., Tk — Tk.o 
will play an important role in constructing the threshold 
equation. We now consider relative positions of the ma- 
jor peaks. For the field \g\iJ,BB = A^, the major peak 
of Pd^ioj) appears at 2A^ above the Fermi level, while 
that of Pdiioj) is located at the Fermi level. In this case, 
no separation will be observed because the major peak 
of Pdi(w) for a reversed bias has the same position with 
that of the forward bias. The spectral densities for for- 
ward and reverse bias are mirror-reflected with respect 
to the vertical axis. 

We assume that the separation may be clearly ob- 
served when the major peak of pdi{uj) appears at a 
position lower than A^ below the Fermi level, i.e., 
\g\pBB + A'^{TK) > 3A^(rK,o)- If we combine this with 
the effect of Kondo temperature mentioned above, we 
write the threshold equation as \g\pBBc + A^(rKm) = 
3A;^ s[l + (rKm-rK,o)/4rK,o], where 

^K,s — A^(rK,o) 

and Tkih is the maximum Kondo temperature at thresh- 
old. At Tk = Tk,o, Be 2A^ g/\g\pB, which corre- 
sponds to Be = 2.4T for A^ g = ll^eV and \g\ = 
0.16[3]. The coefficient 3AjJ; 5/4Tk.o has been introduced 
phenomenologically. Remarkable agreement with the ex- 
periment is shown in Fig. 3, if we use the experimental 
value Tk^o = 0.25K[|. 

We successfully explained the experiments using Eq. 
(1). Now, it is right time to derive Eq. (1) by calculating 
the retarded Green's function Gj^|(w) under nonequilib- 
rium conditions. We use the NDT that provides a new 
technique to calculate G (w) [1| , which is expressed as 

^drftrM = (cdcr|(t^ + irj - L)-^\cda) in the Heisenberg 
picture 6], where L is the Liouville operator defined by 
hA — HA — AH where H is the Hamiltonian, 77 is a 
positive infinitesimal, and the inner product is defined as 
{A\B) = {{A,B'^}), where A and B are the operators of 
the Liouville space, the curly brackets denote the anti- 
commutator, and the last angular brackets represent a 
nonequilibrium average Q. If we define the elements of 




FIG. 3: Relationship between Tkiu and Be at the splitting 
threshold. 
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the matrix M as My = (ej|zI+iL|ei), where z 



-iu+r] 



and the operator Bj is one of the bases spanning the Liou- 
ville space, the retarded Green's function is represented 
by iGjdai^) = {cda\M-^\cda) = (adj M)<id[det M]"!, 
where (adj M)y denotes the cofactor of the jz-element 
in the determinant of M. 

The NDT is initiated by constructing a complete set 
of dynamical bases {ej\j = 1, 2, • • • , oo} spanning the re- 
duced Liouville space in which the dynamics of the op- 
erator Cda{t) effectively describes the Kondo processes. 
We simply extend the bases for the single-impurity An- 
derson model used in the previous study[5| to the case 
of two metallic reservoirs described by the Hamilto- 

manH = T,k,a.a=L,R'^k4lc'^^ + J2k,a.a=L,Ri'^kdcl^c'^„ + 
^kdCtlcda) + Ea=±l(e<i - <^9\fJ'B\B)cl^Cda + Und^Tldl, 

in which L and R denote the left and right leads, re- 
spectively, and the subscript k denotes the quantum 
state of the metallic leads and ridi ~ c^j^^Cdi- Then, 
the basis operators comprise the following five sets: (a) 
Sj^ = {c^-^\k = 1,2, •••,00} for describing the move- 
ments in the left noninteracting metallic lead, (b) = 
{(^nrfjC^I j/c = 1,2, •••,00} for describing the annihila- 
tion of a spin-up electron in the left lead combined 
with the density fluctuations in the spin-down electron 
at the dot, (c) and (d) the sets similar to (a) and (b) 
for the right lead, i.e., Sj^ and S^, respectively, and 

(e) Sd = {5jdi^Cdi,Sj^['cdT,CdT,Sjj^^CdT,jdi^Cd^} for de- 
scribing the dynamical Kondo processes at the dot. 

If we symmetrically arrange the bases such as Sj^, S^, 
Sd, S^, and to construct the matrix M, the following 
matrix of nine blocks is obtained: 
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where the blocks Mc(, M-dL and M^^, and 'M.Ld and M^d 
are 5 x 5, 5 x 00, and 00 x 5 matrices, respectively. Blocks 
M^L and M/j/j are the 00 x 00 matrices that are con- 
structed by the sets of bases describing the left and right 
leads, respectively. Since no direct coupling exists be- 
tween the left and right leads, zero matrices occur at 
the two corners. The structure of each block is similar 
to that of the corresponding block of the matrix for the 
Anderson model considered in the previous studyQ. 

The infinite-dimensional matrix can be reduced to a 
finite-dimensional matrix using the Lowdin's partition- 
ing technique 3, Q. It is performed by solving the eigen- 
value equation for the matrix M^^^, such as MideC — 0, 
where C and are the infinite-dimensional column vec- 
tors. The column vector C is partitioned into three 
parts, i.e., = (C^C^C^), where T, L, d, and R 
denote the transpose, left lead, dot, and right lead, re- 
spectively. Then, the equation for is obtained as 
(Md - MLdMllMdL - MRdM^RMdR)^ = MdC = 
0. The reduced 5x5 matrix M^ contains information 
on the many-body dynamics of the spin-up electron. Ob- 
taining M.d is practically possible if the 00 x cx) matrices 



Mll and M/j^j are block diagonal^. The second and 
third terms appear as self-energies in M^ after reduction. 
The matrix M^ is expressed as 
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g\^B\B. All the matrix 
additional self- 

iTj^{ijj)), where 



where uj = 

elements, other than Uj± and U*j±, have additional self 
energy functions iYjij{u!) = (3ij{iYj^{ijj) 
So(w) is the self-energy for the Anderson model aX, U = 
Opilj. The coefficient /333 = 1 and others, Pij — (3ji, 
are given by /322 = 024 = Pa > /3i2 = > Pii = 
> /3i5 > 1^25 = /945[ll|- In equilibrium at half-filling, 
however, /3y = 1/4 except Ps^. We use that jll = Irr 
and 7j- = 7j+ . The latter will be discussed below. 

In particular, it is noteworthy that the zeros of the de- 
terminant of Mc( are w = 0,±[y^j^ + {'-fLfi—jj± )'^]^^'^, and 
zLU/2 in the large- C/ and atomic limit. The second zeros 
become ±7^^ under appropriate amount of bias. This 
implies that there are additional resonant tunneling levels 
at ±7ii if we treat the system with two separate leads. 
One of the effects of the magnetic field is the Zeeman shift 
appearing in tD in the diagonal elements. We use lu to rep- 
resent uj — ed — U {udi) hereafter. The Zeeman shift results 
in a nonvanishing value of 1/2 — {udi) that induces asym- 
metry in the spectral density by the nonvanishing imagi- 
nary part of the elements Uj± , which can be expressed as 
Uj^ = (C//4)[l + z2(l-2(n,i))0-±)/((<5j±)2)V2]. The 

asymmetry maximizes the new resonant peak appearing 
at = g^isB + 7ll, as shown in Fig. 2. Therefore, 
the field-independent part of the splitting g is merely 
jll- The real part of Uj± determines the positions of 
the J7-peaks. The coefficient 1/4 was determined from 
the analysis performed at the atomic limit. It is l/\/2 
times smaller than that of the single impurity Anderson 
model[l|. 

Since 7LL = + ykd^c^MlijJ^'^ ^ j+l^]) , 

its direct calculation is difficult [l^. We have tried direct 
calculation of 7ll in our previous workp^. We have 
obtained that jll oc ln{D/Zr). However, we failed to 
derive a correct prefactor of lii(D/ZT) in the previous 
work, which is necessary to explain the experimental re- 
sults quantitatively. Even though the NDT provides us 
all qualitative features of the Kondo peak splitting phe- 
nomenon, we need correct prefactor of lii(D/ZT) for the 
quantitative comparison with experiment. An indirect 
way employing the exact result by Bethe Ansatzfl(?| make 
it possible to derive an appropriate expression of jll- 

The wavefunction renormalization Z can also be ob- 
tained by calculating the spectral weight of pd-\{uj) at 
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FIG. 4: Coherent transport through the new resonant tun- 
neling level when a spin-up electron hops into the dot. 



Lu — 0. Under equilibrium with particle-hole symmetry, 
it is given by Zs = [1 + U^/MiIl + iIr)}]~\ where 

Kondo regime, Zs gets separated into two parts, i.e., 
= A-/Il/U^ and = A-flJU^. The former is 

the spectral weight remaining at the Fermi level, which 
does not contribute to the Kondo peak splitting, and the 
latter, which is responsible for the Kondo peak splitting, 
transfers to the new resonant tunneling level at lo — 7^,^ 

and gives rise to the additional shift ^ Zg^U/2 to the 
one by Zeeman effect. 

As one can see from the operator expression of jll, 
it is independent of the Coulomb repulsion. There- 
fore, U appearing in the expression of 7^^ should 
come from the averaging process and may be expressed 
by Z. The theory using Bethe ansatz yields Zs = 
(A/TT)yJU /rexp[— 7r[//4r] for the symmetric Anderson 
model in the Kondo regime[l^. Then, U is written as 
U = (4r/7r)ln(i?/Zsr), where D = Zsrexp[7r;7/4r]. 
For the Anderson model with two separate reservoirs, 
however, the effect of the two reservoirs changes U 
and F into U/y/2 and 2T, respectively. Therefore, 

Zs = (4/7r)y'c//2V2Fexp[-7r[//8\/2F] for the SET 
in the Kondo regime. By using the expression for 
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U from the Zs for the SET and removing the sub- 
script S to take the asymmetry into account, we ob- 
tain 7LL = {^^/2T/^T)^Jz^h).{D/ZV), where D = 

ZsFexp[7rL//8\/2F]. This is just Eq. (1). 

It is interesting and meaningful to observe the mech- 
anism of electron transport through the new resonant 
tunneling level in an SET. This mechanism appears in 
the operator expressions of ^lr- The motion of the spin- 
down electron is governed by the operator jjj^], as 
described in Fig. 1 (a) , while spin- up electron moves into 
the dot from the left or the right lead. Figure 4 shows 
the process making current when the chemical potential 
of the right lead is higher than that of left. Since the 
movement that costs energy U is not allowed, a possible 
operation that constitutes a current is the one in which 
a spin-up electron moves from the right lead to the dot, 
makes a singlet coupling with a spin-down electron on 
the right lead, and moves together toward the left lead. 
Then, the spin-down electron changes its partner that is 
a spin- up electron on the right lead to make a singlet. 
This new singlet moves to the left and the process is re- 
peated. It is noteworthy that the singlet state is retained 
during transport and the partner of singlet is changed, 
similar to that in superconductivity. 

The motion of the spin-down electron of the 
operator in 7^^ = (Efc^l^l^^c^t + 

n:fc« )4Tbd?.jJi'']) is Fig. 1 (a) with - sign, which 
corresponds to the net flow of the spin-down electrons. 
When the bias increases, 7jt approaches ^lr and the 
spectral weight of the new resonant peak also increases. 
However, at equilibrium, 7jt vanishes and the new res- 
onant peak disappears. 

In conclusion, we report the existence of a new res- 
onant tunneling level in an SET with a strong correla- 
tion. This new resonant tunneling level is activated only 
under nonequilibrium conditions, and it explains all the 
features of the Kondo peak splitting phenomenon mea- 
sured by Amasha et al.[i\. Further, we show in Fig. 4 
the coherent transport mechanism through the dot using 
the new resonant tunneling level. 
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